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Abstract 

We solve a new chiral Random Two-Matrix Theory by means of biorthogonal polynomials for 
any matrix size N. By deriving the relevant kernels we find explicit formulas for all {n, fc)-point 
spectral (mixed or unmixed) correlation functions. In the microscopic limit we find the correspond- 
ing scaling functions, and thus derive all spectral correlators in this limit as well. We extend these 
results to the ordinary (non-chiral) ensembles, and also there provide explicit solutions for any finite 
size N, and in the microscopic scaling limit. Our results give the general analytical expressions 
for the microscopic correlation functions of the Dirac operator eigenvalues in theories with imag- 
inary baryon and isospin chemical potential, and can be used to extract the tree-level pion decay 
constant from lattice gauge theory configurations. We find exact agreement with previous compu- 
tations based on the low-energy effective field theory in the two special cases where comparisons 
are possible. 



1 Introduction 



In this paper we propose a chiral Random Matrix Theory of two coupled hermitian matrices, and solve 
the theory exactly at finite matrix size A^. In the N ^ oo limit we find the associated microscopic 
scaling regime, and derive closed expressions for all spectral correlation functions. It has been shown 
recently [li |2j that these spectral correlation functions describe correlations between the lowest-lying 
Dirac operator eigenvalues in the /? = 2 Dyson universality class of spontaneous chiral symmetry 
breaking [3], once subjected to an external vector field. The effect of the external vector source on the 
chiral condensate was considered in ref. [1] and very recently spectral sum rules have been derived 
directly from the effective partition functions [5] . The chiral Random Two-Matrix Theory generalizes 
the original one-matrix problem of chiral Random Matrix Theory [HI [71 El E] (the microscopic hard- 
edge problem with a certain determinant measure) and is also related to the chiral non-Hermitian 
Random Two-Matrix Theory [TO] for QCD with a chemical potential by setting the chemical potential 
to be purely imaginary. 

The solution of the chiral Random Two-Matrix Theory is of interest in its own right, but it is 
also of importance in the context of lattice gauge theory. This is due to the fact that the spectral 
correlation functions of the chiral two-matrix problem can be used to measure the pion decay constant 

in a manner analogous to how the chiral condensate S can be extracted by means of the original 
microscopic solution of the chiral Random Matrix Theory of just one chiral matrix [11] . The idea is to 
compare the analytic finite- volume scaling functions of Dirac operator eigenvalues to lattice simulation 
data. Fits to these data yield accurate determinations of the two low-energy constants of QCD, S and 
F^. 

Before we describe the relevant chiral Random Two-Matrix Theory, we will briefly outline how this 
two-matrix problem arises in the context of particle physics and how F.,^ enters the game. The coupled 
chiral two-matrix problem turns out to describe the microscopic limit of QCD Dirac operator spectra 
close to the origin. The inclusion of one additional matrix in the way we shall describe below allows us 
to determine the correlations between Dirac operators with different imaginary chemical potentials. 
In detail, we shall consider two different Dirac operators, 

Z)i4") ^ [0(A)+i/xi7o]V'l"^ = iA^'^Vl'^^ 

Z)2V?) ^ = iA(")4"^ (1.1) 

The interest in studying these Dirac operators has already been hinted at above: The correlation 
functions between eigenvalues of Dirac operators at different imaginary chemical potential is extremely 
sensitive to the value of Fj^. In the so-called e-regime of QCD this dependence on is uniquely 
determined by the pattern of chiral symmetry breaking. As we will show below the microscopic 
correlations between the two different sets of eigenvalues Ai and A2 only depends on the difference 
— A*!- can then for convenience set /i2 = —^1 = ^ and study the theory with an imaginary 
isospin chemical potential /i. Alternatively we will also consider the case where /ii = and 112 7^ 0. 
Then Di is the usual Dirac operator of a given gauge configuration A^{x) and is what is typically 
used in lattice QCD simulations. This allows one to extract the low energy constants from existing 
simulations. 

A powerful way to extract information about the lowest-lying Dirac operator eigenvalues is through 
the effective low-energy theory of Goldstone bosons, the associated chiral Lagrangian [12]. In the 
e-regime this effective low-energy field theory receives, to leading order, only contributions from mo- 
mentum zero modes of the Goldstone bosons. In addition, the topological index z/ of the given gauge 
field configuration now plays a crucial role, and it is advantageous to consider the partition function in 
a sector of fixed topological charge v. For a theory of Nf quark fiavors, we let B denote the A'^^ x Nj 
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diagonal matrix of imaginary chemical potentials, i.e. for the two Dirac operators of eq. (jl.ll) a diag- 
onal matrix of /ii's and/or In the effective SU(A^/) low-energy Lagrangian the fundamental field 
is taken as U{x) £ SU(A'^^). The translation table from quark couplings to couplings in the effective 
theory is well known. The quark chemical potentials is of vector-kind, and only in the time direction. 
The coupling to external imaginary chemical potentials described by the matrix B is then achieved 
by replacing the time-derivative, 

doU{x) VoU{x) = doU{x) -i[B,U{x)]. (1.2) 

In both the more conventional "p-regime" and the e-regime the leading-order effective Lagrangian is 



-Tr 
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VoU{x)VoU^x) + diU{x)diU\x^ - -Tr \MU^x) + M^U{x^ , (1.3) 



where A4 is the (diagonal) quark mass matrix of size Nf x Nf. This effective theory simplifies even 
further in the e-regime, where to leading order only momentum zero modes play a role. At fixed gauge 
field topology the partition function then reads 

Z,= [ dU (det [/)-eiVF.^Tr[t/,B][C/t,B] + isyTr(A^tf;+^,t;t) ^ 
Ju{Nf) 

where the integration is performed over the full unitary group. Note that -F^r enters in the dimensionless 
combination VF^/j^'j. In the case where all chemical potentials are equal (imaginary baryon chemical 
potential), B is proportional to the identity and all and i^7r-dependence trivially drops out of the 
effective low-energy Lagrangian as expected since the Goldstone bosons do not carry baryon charge. 
However, when different chemical potentials enters there is a unique dependence on F^. The correlation 
function between the spectra of Di and D2 have a dramatic dependence on F^. In [H [2] the two 
point function was obtained for zero and for two flavors directly from the low-energy effective chiral 
Lagrangian. Here we generalize these results using a chiral Random Matrix Theory. The chiral 
Random Matrix Theory approach to the e-regime in this case is a two-matrix problem. We need two 
matrices because we wish to describe two different kinds of operators, Di and D2- With this extension 
of the original chiral Random Matrix Theory we will derive all (n, fc)-point correlation functions of 
the operators Di and D2- The advantage of the chiral Random Two-Matrix Theory approach over 
the direct calculation using the low-energy effective chiral Lagrangian lies in the ease with which 
the calculations generalize. We will verify explicitly that the chiral Random Two-Matrix Theory 
reproduces, in two special cases, the results obtained previously through the low-energy effective 
chiral Lagrangian [H [2] . 

In this paper we will also find the corresponding solution of the non-chiral Random Matrix Theory 
with two coupled matrices generalizing the results of Eynard and Mehta ^13j to include determinants 
(matrix analogues of Dirac determinants). This is of relevance for Dirac operator spectra of the 
analogous 3-dimensional gauge theory [14^ [T5] . It may have applications in condensed matter physics 
as well fW\- 



We round off this introduction with a few remarks on the case of real chemical potential. The 
spectral properties of the Dirac operator change drastically if we were to consider instead ordinary 
real baryon chemical potential. While both Dirac operators (jl.ip are anti-hermitian (and all A^"^ thus 
real), this is not the case for real chemical potential. The eigenvalues in that case lie in the complex 
plane. The precise manner in which the eigenvalues spread into the complex plane when a real baryon 
chemical potential is turned on has been computed in the e-regime |10j by a random matrix approach 
and subsequently verified by a low-energy effective calculation [17] . 
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A Random One-Matrix Theory describing real chemical potential in the low-energy limit was first 
suggested by Stephanov [18] . As shown in p[Oj , the way to find an appropriate eigenvalue representation 
is to start out with a two-matrix extension, of which one matrix is anti-hermitian, and the other 
is hermitian. The eigenvalues of the combination are complex, and their microscopic distributions 
coincide with the low-energy Dirac operator eigenvalues at real baryon chemical potential. Much work 
has recently gone into the understanding of this low-lying Dirac spectrum [TUl [201 ED Eni [23 El 123] 
(for numerical evidence from lattice gauge theory simulations see, e.g., ref. [23j)- As stated earlier the 
Random Matrix Theory at real baryon chemical potential introduced in [lOj is related to the chiral 
Random Two-Matrix Theory for imaginary chemical potential we propose in this paper. We will point 
out the similarities and differences below. 

Our paper is organized as follows. In the next section we introduce the chiral Random Two-Matrix 
Theory, and proceed to derive an eigenvalue representation. This leads us to the use of biorthogonal 
polynomials, which we identify explicitly for finite size of the matrices in section [3l There we present 
our results for all spectral correlation functions in both the quenched case, the partially quenched 
case, and the full case (appropriate definitions of these terms will be given below). We also derive 
the microscopic scaling limit in all three cases. In section [J] we turn to the non-chiral variant of this 
theory, and find also there the explicit solution both at finite and in the analogous scaling limit. 
Section [5] contains our conclusions. Several technical details are collected in 3 Appendices. 



2 The Chiral Two-Matrix Problem 

The chiral Two-Matrix Theory with imaginary chemical potential is described by the partition function 

Nf 

zi^^'\{m}) ~ / d^d"^ exp [-A^Tr (^>t$ + ^rt^^j j| det[V{nf) + nif] (2.1) 

/=i 

where the anti-Hermitian operator D in the determinant is given by 

Both $ and ^ are complex rectangular matrices of size A^ x (A^ -|- v), where A^ and u are integers. As 
shown in ()2.ip we restrict ourselves to Gaussian weight functions for both ^ and ^ in this paper. The 
partition function [10] at real chemical potential is obtained by taking all fij purely imaginary. Because 
this changes the Hermiticity properties of the Dirac operator the eigenvalue correlation functions at 
real and imaginary chemical potential are not related by a simple analytic continuation /ij i^f. 

We stress already here that ^ f and m f are (dimensionless) quantities that a priori have no relation 
to either chemical potential or quark masses in QCD. However, in the large-A^ scaling limit that we 
shall discuss below they do translate into the familiar low-energy parameters of QCD by means of 
/ij = — > ^'jFI^V and rh = 2mN — > mT,V. With this in mind we refer to Nf as the number of 
flavors (here we consider only Nf > 0). 

We will consider the eigenvalues of the two matrices Pi = T>{^i) and P2 = ^(/f^2) for arbitrary 
(real) values of //i and ^2- The case of imaginary isospin chemical potential can be obtained from this 
by setting ^2 = — /^i = A*- To obtain the eigenvalues we flrst change integration variables to 

$1 = $ + 

$2 = ^ + /i2^- (2.3) 
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We will restrict ourselves to the case where all determinants in (12. ip contain T>i or I?2- Their numbers 
are denoted by A^i and N2, with masses rufi and mf2, respectively: 



^^({mi}, {ma}) ~ / d<^d^ H det[V, + m/i] H det[P2 + ^^2] e-^Tr(4.t$+vi>t.,) ^^.4) 

/1=1 /2=1 

with Nf = Ni + N2 ■ After the change of variables (j2.3p the partition function becomes 

Zl'^'^imi}, {m2}) - d^id^2 n det[Pi + mji] [] det[p2 + m/2] e-^TrV{*i,*2) (2.5) 

/1=1 /2=1 

where each of the two operators are of the standard chiral form, 

/ \ , , 

^1 = I .^t n I (2-6) 







and 



However the two matrices and $2 are now coupled: 

$2) = Ci«>i$i + C2$^«>2 - d ($1$2 + (2.8) 

with the abbreviations 

Cl = {l + fil)/6^ 

C2 = (l + ^?)/5' 

d = (1 + ^1^2)/'^^ 

S = /i2 - m , (2.9) 

used in the following. This is the chiral analogue of the Two-Matrix problem introduced in the classic 
papers of refs. [251l26i| . We solve this chiral theory by the biorthogonal polynomial method, analogous 
to what is described in ref. [26]. At this stage, after changing to variables $1^2 we could in principle 
allow for the following more general class of weight functions o 

V{^i,^2) ^vi(<^\^i^ +V2[^l^2) -d(^^\^2 + '^l<^i) , (2.10) 

where vi^2 are polynomials. It remains an interesting open challenge to establish universality of the 
present results in the large- limit. A generalization of the universality proof in ^ i9j is non-trivial 
due to the different structure of recursion relations for biorthogonal polynomials in the two-matrix 
theory. 

The first step now consists in finding an eigenvalue representation for the partition function (j2.5p . 
To this end we follow the standard procedure [6j and write 

$1 = UiXV^ 

^2 = U2YV^ , (2.11) 



^Switching to higher order polynomials before the change of variables leads to higher order terms in <E>|<1>2 and its 
conjugate. 
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where C/1,2 and Vi^2 are unitary matrices (with Vi^2 restricted to U{N + i/)/[U{l)]^~^'^ [7]), and X and 
Y are diagonal matrices with real and non- negative entries Xi respectively yi. Under this change of 
variables the measure (including the determinants) becomes 

A^i N / Ni \ 

d<^i n det[Pi + m/i] = dUidVi J] dxjxj" [] i^l + ^N{{x^}f 

fi=i i=i V /i=i / 

d^2 n det[P2 + m/2] = dU2dV2 n dyfyf'' Y[ (yf + m^) Aiv({y'})' (2.12) 

/2=1 i=l \ /2=1 / 

where we have left out the factors of rrij coming from the zero modes of the determinants. Here 
and Vi^2 are integrated over the Haar measure, and A^v is the Vandermonde determinant. Because of 
the coupling between $1 and $2 the integration over Ui^2 and Vi^2 is not triviaH. Using left and right 
invariance of the Haar measure the required integral reduces to a form that has been evaluated in ref. 

m- 

[dUdVexp[dNTriVXUY) + h.c. ] = ^'a'^?/ n^A ^1/ 2., ^ (2-13) 

up to an irrelevant normalization factor. Here Iu{x) is a modified Bessel function. 
Putting everything together we thus finally arrive at the eigenvalue representation 

... . roo N / m N2 

Zr\{mi},{m2}) = Af{{mi,2}) U \ ^^^dy^ixmY^' U i^' + ^h) U iv' + ^h) 

i=l \ fl=l /2=1 

X AN{{x^})Ajv{{y^}) det [U2dNx,yj)] e'^Ef (2.14) 

where the factors of m'j from the zero modes have now been absorbed in the normalization constant J\f. 
In contrast to the one-matrix theory we cannot extend the integration range to —00, as the integrand 
is no longer even in Xi and yi. The resulting expression (j2.14p has a quite striking resemblance to the 
non-hermitian theory describing the Dirac operator spectrum in the presence of real baryon chemical 
potential [10]. In that case the explicit Bessel function in the measure is of the kind K^, rather than 
the of eq. (|2.14p . In the eigenvalue representation (j2.14|) we can apply the technique of biorthogonal 

polynomials. For that purpose we seek two sets of monic polynomials P^^\x'^) and Qn^^\y'^) that 
are biorthogonal, 

~^dxdy w^''f\x,y) Pr^')(x2) Q^^'\y^) = h^n'hnu , (2.15) 


with respect to the weight function 

A^i N2 

w^^f \x,y) = {xyy+^ Hix'^ + m}^) {y"^ + m^) lA'^dNxy) 6-^^'^'''+''^^'^ (2.16) 

/1=1 /2=1 

Note that the order of arguments in w^^f \x, y) is important, as x only appears in the A^i mass terms, 
while y only appears in the A'^2 mass terms. Only in the special case of 112 = — A*i) -^1 = -^2 and masses 
pairwise matched, mji = m/2) is the measure symmetric in x and y. This is in particular true for the 
quenched case with an imaginary isospin chemical potential. 
After the usual rewriting of the Vandermonde determinants, 

A^({x2}) = detlxf'-'^] = det\p^^(\x'f] 



■^Of course this is only true for /ii,2 7^ 0. 
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A^({y2}) = det[yf^-^)] = det [qfjiyf)] (2.17) 

in terms of polynomials in squared variables one can use the biorthogonality condition ()2.15p to 

perform the integral (|2.14|) in terms of the norms hl^^^ giving zl^^^ = NIJJ^^q hil^^\ The result 

can be expressed in terms of the kernel Kn{x, y) of the quenched polynomials Plf\x'^) and Q^^\y'^), 
defined in eq. (j2.20p below. Details of the derivation following pO] are given in appendixO leading to 

ATi N2 jrN+Ni-l ^(0) 

/lii /iii ^N^[{mi})AN2{{mi}) 



det [ 



for A''i > N2. Here we have specified the normalization leading to N{{mi^2]) = 11/1=1 "^/i n/§=i "^/2- 
We are not only interested in the partition function itself, but mostly in the spectral correlation 
functions. The biorthogonal polynomial technique has been adopted to this purpose in refs. [28 ^1131129] . 
We review it briefly here, and in particular we make use of the very compact and general formulation 
of |13j . The measure considered in that paper is of the ordinary unitary theory which through the 
Itzykson-Zuber integral [25] gives rise to an exponential coupling between the two sets of eigenvalues. 
In our case this is replaced by the coupling through the Bessel function Jjy as in eq. (12.150 . However one 
can straightforwardly follow through the same steps as in flS] with this replacement. In general four 
different kernels are required. Three of the kernels also depend on the generalized Bessel transforms 

X^u'\y) - rdxu;W)(x,y)Pr^)(x2) 

JO 

xf^'^(x) - rdyw^''f\x,y)Qf^\y'). (2.19) 



Note that in general the two functions are not polynomials. In particular, they are proportional to 
the part of the weight which is not integrated. 

The four required kernels in terms of the biorthogonal polynomials and their transforms are: 

fe=0 '^^ 

<^)(xi,X2) = E ' , = / dyw('^fHx,,y)K^^''^\y,X2) (2.21) 

k=0 Jo 

H'^^yuy.) = j: ^'"''^y\^'^''^y^^ = fdx.(-.)(x,..)i^f^)(.i,.) (2.22) 



Mjf^-^(x,y) = E ^t^' = r dudvw^''f\x,u)w^''f\v,y)K^^'\u,v).{2.22,) 



k=o 

Note that also here the order of arguments is important. We have chosen the notation and yi to 
indicate the nature of the arguments so that Xi is an eigenvalue of associated with Pi, while yi is an 
eigenvalue associated with 2?2- 

All spectral correlation functions can be derived from the four kernels above. In the special case 
of a symmetric weight, w^^f\x,y) = w^^f\y^x)^ the polynomials coincide, Pjf^^^ = Q^^^\ leading to 
coinciding transforms x~^^^ = X^^'^ ^"^^ thus to a matching of the kernels H^^\x,y) = H^^^\y,x). 
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We call the correlation function of a set containing n eigenvalues of Di and k eigenvalues of P2 an 

(n, k)- correlation function 



iV!2 



1 



{N-n)l{N-k)l2l''f\{m,},{m2}) 



(2.24) 



N 



N 



n dx, n dyjdet[w^''f\x„yj)\AN{{x^})AN{{y^}). 

i=n+l j=k+l 

Note for later comparison that these include both the connected and disconnected parts. 

Eynard and Mehta [T3] have proven a generalized "Dyson Theorem" for all correlation functions 
of random multi-matrix theories coupled in a chain. Adapted to our chiral Two-Matrix Theory it can 
be stated as follows: 



<i)({^}n>{y}'^) 



det 



K 



(Nf) 
N 



{Ujl 5 -^12 



H'N'\yji,yj2) 



The simplest example, the spectral density itself equals 



R 



(Nf) 
(1,0) 



{x) = Hi^^\x,x) 



(2.25) 



(2.26) 



and similar for i?|^^j(y). In the symmetric case (isospin) the densities are the same for the two 
operators 'D{fi) and T>{—fj,) as is obvious from the symmetry ^ in the partition function (j2.ip . 

More examples will be given in the next section. If all eigenvalues of one kind are integrated out 
this reduces to the known result for the one- matrix theory (up to a /i dependent rescaling), as can be 
seen by generalizing the argument from [28] for the unquenched non-chiral theory: 



i?l^{L,({x}„) 



det 



H 



(Nf) 
N 



{X'i, Xj) 



(2.27) 



and similarly for -R(^io^fc)({y}fc)- 



3 Explicit Computations 

In this section we explicitly determine the biorthogonal polynomials fulfilling the condition (j2.15p . It 
is instructive to consider first the quenched case Nf = where there are no fiavor determinants in the 
integration measure. This is not only because matters simplify then, but because the general solution 
can be expressed in terms of quenched quantities. Thus after having gone through the quenched case 
in detail, we present the general solution for any Nf. We also take the large- scaling limit, and thus 
make contact to the field theory calculations of refs. [H [2]. 



3.1 The quenched case 

The analog of quenching in the Random Matrix Theory formulation is easily achieved: We simply 
remove the explicit determinant factors in the measure. The theory remains chiral, and we still 
perform the change of variables ()2.3p . The weight function (j2.16p then simplifies considerably and the 
orthogonal polynomials can be readily obtained as shown in Appendices El and |Bl The result is 

p„(2;2) = nli-NTCiyL^^^NTCix^) 

Qn{y^) = n\{-NTC2)-"Ll,-\NTC2y^) (3.1) 
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with 

T = 1 - <f/{ciC2) (3.2) 

from Appendix [B] and the normahzation constant is given in (jB.lip . Here and in the following we 
drop the superscript {Nf = 0) for the quenched quantities. The transforms (|2.19p can also be found 
from ([Ri]) along with pIOl) and ([RTO]) 



These are proportional to the polynomials in eq. (j3.ip . 

Having the required functions, we can now write down the related kernels. These can be obtained 
from (|B.5P which gives 

KN{x,y) = 1. (fc^^), [-^) Ll> [Ntcix jL',' [NTC2y ), (3.4) 

H^{x,,X2) = 2(iVrci)'^+ixr+ie-^-i-? ^ (fcT^^^'^ (iVrcix?) {Nrc,x: 

N-l 



k=0 



HN{yi,y2) = 2{Nrc2rV2-'^'e-^-^^y'^ Yl JT^A^ (^^^^y?) ^^2^1 

fc=0 l'*^ + '^J- 

M^ix,y) = riNrdYixyr+'e-^^^'^^^'^'^^y'^ J: 4^^) (iVrcix^) 4^) (iVW) 



fc=0 



We now take the microscopic scaling limit where N ^ oo and fi'j = 2fi'jN and m = 2mN are kept 
fixed. We define microscopic variables by 

X = 2xN 

y = 2yN . (3.5) 
Making use of the asymptotic behavior of Laguerre polynomials, 

lim 4")(x2) ~ n"(nx2)-"/2ja(2V^x) (3.6) 

n— >oo 

where Jk{x) is a Bessel function, we can turn sums into integrals '^/NJ2k=o ~^ fo with t = k/N 

f^Jkf^X'^y'^'^' (^-^-') (^-^^') -ji^l'dte^^'^/^MxVDMyVt) . (3.7) 

Here s can have the values 0, ±1 and Cx and Cy can be either of ci or C2. We have also used 



with 6 = fi2 — fj-i- In this scaling limit the kernels therefore take the form 

lim KN(x = —,y = -^] = 2^''+^N^''+\xy)-^6-^I+{x,y) (3.9) 
Af-.oo 2iV 2N J ^ ' ^ ^ ^ ' 
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lim Hn I X 



lim Mm x 



X _ y \ 
2N'^ ~ 2NJ 



X y \ 

— :,y = — ' 



2Nx''+^y-''l\x,y) 
lim Hn (y = X 

N^oo \ 



X \ 

2n) 



2N' 



2N J 



where we introduced the abbreviated notation 

1 /•! 

2 
1 

2 

The weight becomes 



2N' 

^--\xyy+H^l~{x,y) 



(3.10) 



(3.11) 



J^{x,y) 
I^{x,y) 



dt Ju{x\/t)Jy{y\/t) 



xJu^i{x)Ju{y) - yJu+i{y)Jv{x) 



x^ — y'^ 



w«»'(i,S) = Jim»™(^,^) 



(o)f A = ( \ 



v+l 



exp 



2(52 



xy 
P 



using the abbreviation 



T {x,y) 



P 



exp 



x"^ + y2 
2P 



xy 
P 



I {x,y) . 



(3.12) 



(3.13) 



(3.14) 



Notice that the /u's only appear through their difference. Despite the relation ()B.6P between the weight 
and the kernel M we will derive in appendix [B] our large- results of the two differ. The reason for 
this is that the large- limit is not uniform: while we have rescaled the chemical potentials with N, 
also called "weak non-Hermiticity" in the context of real //, another limit exists where the chemical 
potentials are not scaled, including a different rescaling of the eigenvalues with N. In this limit called 
"strong non-Hermiticity" for real /i, which can also be obtained by letting 6 ^ oo, the large- A'^ results 
of weight and AI will match, as suggested by relation (IB.6|) . For more details on this limit we refer to 

m- 

The correlation functions in the microscopic scaling limit are defined as follows for general Nf. 

1 



lim 



-R 



(Nf) 



{x 



:}n,{y 



y 



jv^oo (2iV)"+fc 2Ar-""^" 2N 

Inserting eqs. (j3.9p - (j3.1ip into the definition we arrive at the following result 



'-^^{xi-^, Xi^) i (xjj,yj2. 



(3.15) 



(3.16) 



P{n,k)i{x}n,{y}k) = Y[xiY[yjX ^det 

after taking out common factors. 

As a special case we note that the (l,0)-correlation function reduces to the known Bessel kernel 
taken at coincident points, yielding the /ij-independent 



JI{x) - Ju+i{x)Ji,-i{x) 



(3.17) 



which is the known quenched expression from the one-matrix theory [7j. The ///-independence of this 
result is particularly easy to understand from the effective Lagrangian, where the corresponding B- 
matrix is proportional to the unit matrix, and thus can be shifted away. Similarly, the (2,0)-correlation 
function of this quenched theory becomes 

2 



P(2,0)(^l>^2) = P(l,0)(^l)P(l,0)(^2) - XiX2 



XiJyj^l{xi)Jy{x2) - X2Ju+l{x2)Ju{xi) 



(3.18) 
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which again is /ij-independent, and agrees with the standard expression for the ^/ = theory [7j. 
One easily sees how this generalizes to the the general (n, 0)-correlator in the quenched theory: the 
usual n-point correlation function of the one- matrix theory [HI [9] is recovered from eq. (j2.25p . taking 
out and canceling the prefactors from eq. (j3.10p from the determinant. 

A particularly important quantity is the mixed (l,l)-correlation function P(i^i)ix,y) which reads 

P{i,i){x,y) = P(i,o)(^)P(o,i)(y) + P(M"(^'y)- (3-19) 
The connected part is given by 

p'{°,'\^x,y) = xyl+{y,x)i-{x,y) . (3.20) 

After making the required substitution (5^ — > AVF^^fi'^ and the identification x = A+Sy, y = A_Sy, 
this agrees exactly with what was previously computed from the effective low-energy QCD Lagrangian 
[1]. In that framework, however, the general correlation functions P(°"^ have not been obtained. 

Another interesting case is the quenched limit when p, is not scaled microscopically, as discussed 
above. In this limit the upper right corner of the determinant in eq. (j3.16p vanished due to the 
matching of the weight and large- kernel M. This leads to a factorisation of all correlation functions 
into a product of determinants over H and H. 



3.2 Two AND MORE FLAVORS 

To find the biorthogonal polynomials it is convenient to use the following general expression in the 
eigenvalue representation 

^f'\^') = {-)H^xr''{det[V, + ix])^.^, 

Ni N2 



n dx^dy^{x^y^r+\xj - x') n i^l + ^ji) n + ^'2) 



Zv i=\ \ /1=1 /2=1 

X Afe({x2})A,({y2})det[/,(2dA^x,y,)]e~^S,-i-?+-2J/? (3.21) 

generalising the observation made in the appendix of [29]. As also noted there a similar relation holds 
for 

Q^^'\v') = (-)'=(%)-^(detP2 + ^y)^(iv,) . (3.22) 

The expectation values are taken with respect to the partition function in eq. ()2.5p her^. The prefactor 

{—)'' {ix)~'^ ensures that pjf^^\x'^) is a monic polynomial of order k, and similarly for Q^^^^ [jj^)- Eq. 
(|3.2ip (eq. (j3.22p ) is recognized as being a partition function of one additional Di-flavor (I?2-flavor) 
with purely imaginary mass ix {iy), 

pW■)(^2^ ^ ( )fc(.^)-'^ ^^^'^^'^'^(i"^li'^^'i"^2}) 

' Zf^)({mi},{m2}) 

indicated explicitly in the superscript. In the large- A scaling limit this gives a particularly useful 
representation [30] since compact analytical expressions can be obtained from refs. [2H I22j . where 



^The partition function used here contains k eigenvalues, whereas the exphcit factor of A'' in the exponent remains 
fixed. 
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the involved finite-volume partition functions were derived in the effective field theory framework. 
Alternatively, we can re-derive those expressions from our two-matrix theory by substituting the 
corresponding large- expressions from the previous subsection into eq. (13.23p . and making use of 
the general expression for partition functions given in appendix [Cl 

Instead of inserting the relation eq. (j3.23p into the definition of the kernel eq. (j2.20p a far more 
useful relation exists, giving the kernel directly in terms of a single partition function: 

K'}^il{y,x) = ^—^^{det[Vi+ix]det[V2 + iy])^iNf) (3.24) 

This relation follows in complete analogy to [20] after going through the same steps as in appendix [Cl 
The highest polynomials in K^J^l (y, x) are of degree N, and in order to match we have to divide by 
their coefficient given by pOl [31] 

ui^f) _ / \N^+N2u det[KN+N2+i{im2,imki) PN+N2+i{{imkif) ...PN+NAi^^kif)] fo ok\ 

"Af — \ n'N+N2 T ■ r 7^ T- '■ VT> 7T- \T\ d 7T- vlyT W-^-^J 

det [A^r4.^r2(^m^2,^mfcl) P]si+N2\{'irnkiY) ■ ■ ■ PN+Ni~i(.{i'mkiY)\ 

for Ni > N2 (for N2 > Ni simply exchange the 1 with the 2 quantities as in appendix [Cj. In the 
large- limit the two determinants will cancel, leaving only the quenched norm. In terms of partition 
functions with two additional fiavors this can be written as 

Compared to the polynomial we have added a second fiavor with chemical potential /i2, related to 
the first argument of the kernel. If on the contrary we had added a second flavor of the same kind, 
(det[Pi det [Pi -|- 4x2] ) (Nf) , we would have obtained a determinant of polynomials instead of the 

kernel (see appendix [C]). Only when ^1 = /U2 = the two are related through the Christoffel-Darboux 
identity. 

The functions Xk^^\y) ^^^^ Xk^''^\x) can then be computed from their definition (|2.19p 

xf^'^(y) = r^^^^^^Hx,i/)(-)n^x)- ^-'''^'''''^(^"^^^'^"'^"^^» 



Xk 



^r^^({mi},{m2}) 
Zu {{mi},{m2}) 



Similarly the remaining kernels can be obtained from eq. ()3.26p by transforming one or both arguments 

M^f)t \ /"""j (Nf)f N/ Zi^^^^'^''^^\{mi},ix2,{m2},iy) , . 

H^l^{xi,X2) = / dyw"^ f>{xi,y){-X2y) (^y.) (na (^•2^) 

•^0 hj^^ Zu ({mi},{m2}) 

H^'^fhy, y,) - rdxtf^^^Vx vM ^y^~u 2u'''^''''"^'\{mi},ix,{m2},iyi) 

H,,M,y2) - I dxv. ix,y2)i xy,) ,m ^iN,, ^^^^^^ ^^^^^ (^.29) 

Ml,rA^,y) = r dudv.^-f\x,uW-Hv,y)i-vu)-?^^ 

■^0 h^^ Zu ({mi},{m2}) 

(3.30) 
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Inserting them into eq. (I2.25P this provides the general (n, fc)-correlation function at finite- iV of the 
theory with arbitrary Nf. Since they are expressed in terms of partition functions the microscopic 
scaling limit can easily be taken, as we have mentioned above. 

We now specialize to the physically most interesting case of two flavors A^i = = !> with quark 
masses ruu and m^, respectively. We first focus on the (l,l)-correlation function '^"""(x, y) in 

order to compare with the known result [2j. The first building block we need, the unquenched kernel, 
is given by 



K 



(1+1) 

N+l 



hN 



+1 



■det 



KN+i{imd,imu) KN+i{imd,x) 
KN+i{y,imu) KN+i{y,x) 



(3.31) 

where we used that the quenched kernel is even in its arguments. The second building block is obtained 
from that by integrating twice from eq. (j2.23p 



" dudv (x, u)u'(i+i) {v, y)K^Nti ^) 



^jv+i(x^ + mg)(y^ + mg) 



det 



KN+iiimd,imu) HN+iiimd,y) 
HN+i{x,imu) MN+i{x,y) 



(3.32) 



Taking the integrals inside the determinant leads to the appearance of all different quenched kernels. 
In order to take the large- A'^ limit we can thus use the quenched results eqs. (13. 9p - (|3.1ip . We therefore 
obtain for the two microscopic building blocks 



lim K^^^l\y,x) 



(2Ar)2.+8(^^)-.^-2 / ^ I+{y,i7fiu)I+{imd,x 

(x^ -I- m2)(y2 _|_ ^2-) I 



and 



hm Ml^+l\x,y) 

N — ^■oo 



S^ixyY+^ix"^ + m2)(y2 + rhl) 



'-N+1 

along with the weight 



I (x,y) 



lim ■w^^^^\x,y) 



(2iV)2^'+6 ^ 

(x^ + m^)(y^ + m^) (xyY^^ 



l+{irhd,irhu) J 

l^{x,imu)T^{imd,y) 
I+{imu,iihd) 



■I 



2S2 



(3.33) 



(3.34) 



(3.35) 



(2iV)4 {2N)^^+^' 

Here the norms have canceled in the large- limit, and the explicit mass factors will also drop out 
when multiplying all together in the next step. Using eqs. (j2.25p and (j3.15p we first consider the 
following connected part of the two point function: 



(1+1) conn / * -,\ 

P{i,i) i^'V) 



= lim , 

A^^oo (2A^) 



xy 1^{y,x 



l^Ki^+l\y,x) (^(1+1) (x,y) - Mi^^l\x,y) 
I'^{y,imu)'I~^{imd,x] 



(3.36) 



I+{imd,imu) 



I (x,y) + 



I^{x,imu)T^{irhd,y) 
I+{irhu,imd) 



Substituting 6^ AVF"^^"^ and x = A+SF, y = A_SF, this again agrees with the unquenched result 
from effective low-energy QCD [2J. For completeness we also give the other two building blocks for 
this theory. Following the definition (j2.2ip we have 



H^N+l'{xuX2) 



'dy w^'+'\x,,y)Kl^^l\y,X2 

hN+i{xj + ml) 

^^N^\x2 +rnl)KN+i{imd,imu) 



(3.37) 



det 



KN+iiimd,imu) KN+i{imd,X2) 
HN+i{xi,imu) Hn+i{xi,X2) 
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leading to 



lim H 

N^oo 



(1+1), 
N+l ' 



Xi,X2) 



2Nx'i+\xl + rni) 



For the last block we obtain from eq. (j2.22p 

1 

dxu;(i+i)(x,y2)i^S+!^(yi,x) 

hN+iivl + "T-d) 

^iv"^ ^ (2/1 + rnDKN+i {irud, irUu) 



J°(xi,X2) 



I^jxi, irhu )X+ (imrf , X2 ) 
l+{irhd,irhu) 



. (3.38) 



(3.39) 



det 



KN+i{imd,imu) HN+i{imd,y2) 

KN+i{yi,imu) HN+i{yi,y2) 



leading to 



lim <|)(,„,2) 



yl{yl + ml) 

As a new microscopic result we have the spectral densities 



J+(^i,im„)T°(imrf,7/2) 



. (3.40) 



lim — ^^^+j^(y,y) 
TV^oo 2N ^^^^ 



(3.41) 



given by eqs. (|3.38|) and (|3.40|) at equal arguments divided by 2A^ respectively. Note that the two 
functions coincide. Ai fii = fi2 = they reproduce the known result [9] for two different flavors. For 
a plot of the unquenched eigenvalue density see flgure [TJ 



N^2 m=3 m^=4 |i=5 
Nf=2 ni=3 m^=4 |i=0 




Figure 1: The eigenvalue density for Nf = 2 with imaginary isospin chemical potential (1 = 5/2 = 
and 5 and quark masses m„ = 3, rh^ = 4. The quenched density eq. ()3.17p of the one-matrix theory 
(Nf = 0) is given for comparison. 
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To finish this subsection we give the result for the most general correlation function in the two 
flavor case. If follows by inserting the above building blocks into eq. (|2.25p and taking out common 
factors: 



-n—k 



i=i j=i 

J+ {irhu , irhd)l^ {xi^ , Xi^ ) 



det 
i<iij2<fe 



I+{irhu, imd)'I^{yj^,Xi^) 
-'^^ iVji > ^rhu)I^ (xi2 , irhd) 



(3.42) 



-I+{irhu,irhd)I (yj^ ) Xi2 ) 
-'^^{Vh > irhd)T^{xi^ , irhu) 

I^{irhu, irhd)T^{yj^,yj.,) 
-'^^ iVji > irhu)I°{yj2 , irhd) 



3.3 Partial quenching 

Here we consider a variant of the general A'^^-theory of practical importance. We introduce Nf fermions 
with Dirac operator Di and set ^ui = then seek the spectral correlation functions between eigenvalues 
of Di and D2 for fi2 7^ 0. In lattice gauge theory simulations this is a simple modification. The 
ensemble averages are then taken with respect to a theory which is not coupled to a chemical potential, 
but the observables (expressed in terms of "valence quarks") are. In the Random Matrix Theory 
formulation it can also be introduced straightforwardly. 

The corresponding Random Matrix Theory is based on the partition function eq. (j2.5p with 



N- 



1 = Nf and N2 = and we ask for the mixed correlation function between two valence 

quarks associated with matrices Pi, 2- For simplicity we here restrict ourselves to the physically most 
interesting case of Nj = 2. Here we have two sea quarks with masses rriu and and zero chemical 
potential fii = 0, while the valence quarks have Hi = and /i2 = 5. 

In the previous subsection the results for the most general A'^^ = A^i + A''2 flavors content have been 
given already for the polynomials in eq. (j3.23p and the kernel in eq. (|3.26p from which we can also 
compute the transforms (i239|) . ([MI]), ([2:22]) and 1^2^ . We just have to specify our case Ni = 2, 
A''2 = when applying eq. ()C.5p . We will first focus on p^^^""''^{x, y) in order to pinpoint the difference 
between the 1 + 1 flavor case eq. (I3.36P from the previous subsection, which was also computed in 
an effective field theory approach [2|. Next, we give the general (n, A;)-correlation function for this 
particular flavor combination. 

The first building block for that from eq. (j3.26p is given by 



det 


KN+iiy,imu) 
KN+iiy,imd) 
KN+iiy,x) 


PN+i{-ml) 
PN+i{-ml) 
Pn+i{x^) 


PN+2{-ml) 
PN+2{-mj) 

Pn+2{x^) 


h^^^\x'^ + m1){x'^ + m^) 


det 


" Pn{- 
_ Pn{- 


-ml) Pn+i{- 
-ml) Pn+i{- 


-ml) 

-ml) _ 





(3.43) 



We observe a different structure of determinants compared to the corresponding previous result eq. 
(I3.3ip , marking the difference of partial quenching. The microscopic large- limit is now easily taken 
using the formulas from the previous subsection. The only difference is that we have to Taylor expand 
the last columns in numerator and denominator as they become degenerate at large- A^. The expansion 
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in maps into a Taylor expansion with respect to the arguments and we obtain 



hm K 



(2) 
N+1 



det 

(2iV)2^+8(xy)-^5-2 


I^{y,imu) Jy{irnu) irhuJu+iiirhu) 
I~^{y,imd) Juiirhd) imdJu+iiirhd) 
I^{y,x) Ju{x) xJu+i{x) 


(x^ + m2)(x2 + ?fi^) 


det 


Ju{irhu) irhuJv+i{irhu) 
Juiirhd) irhdJu+iiirhd) 





(3.44) 



The result for the second building block of our correlation function is obtained by integrating eq. 
(f3:i3]l using the definition eq. (12:23]) 



det 



'■N 



I 2 I 2\i 2 I 2\ 

[x +m„)(x +m^) 



HN+i{x,imu) PN+i{-ml) PN+2{-ml) 
HN+i{x,imd) PN+ii-ml) PN+2{-m1) 
MN+i{x,y) XN+i{y) XN+2{y) 



det 



PN{-ml) PN+i{-ml) 



(3.45) 



In addition to the asymptotics already computed we also now need one of the transforms xn- Being 
proportional to the polynomials, eq. (|B.4p . this is an easy task, and we obtain an explicit exponential 
5 dependence from the prefactors. We thus arrive at 



lim M 



(2) 
N+1 



{x,y) 



det 

{xyY^^{Sp' + rn^){x^ + m^) 


I^{x,imu) Jviiniu) irhuJu+i{imu) 
I^{x,imd) Jv{irhd) irhdJu+iiirhd) 
I-{x,y) e-^y^My) e-^'/'GAvJ) _ 


(2Ar)2^+6^2 


det 


Juiirhu) irhuJu+iiirhu) 
Juiirhd) irhdJu+iiirrid) 





where we have introduced the notation 

Gu{y:5) = yJu+i{y) + 5"^ Ju{y) ■ 
Together with the asymptotic of the weight function, 

\x,y) = ^NT^'^ 



xy_ 



we obtain the partially quenched two point function 



(2) conn , ^ ^\ » ^ 

P(i'i) {x,y) = xy 



I+{y,irhu) Ju{irhu) irhuJu+i{irhu) 
I+{y,irhd) Ju{irhd) irhdJu+i{irhd) 
I~^{y,x) Ju{x) xJu+i{x) 



-I^{x,irhu) 
-l^{x,imd) 
i'{x,y) 



Ju{irhu) 
Juiirhd) 

e-''"Ju{y) 



(3.46) 
(3.47) 

(3.48) 



irhuJu+iiirhu) 
irhdJu+i{irhd) 
e-^'/'GuiyJ) 



Ju{irhu) irhuJu+i{irhu) 
Juiirhd) irhdJu+iiirhd) 



(3.49) 

We note that this result is of the structure expected from the replica limit of the Toda lattice equation 

It is ideally suited for extracting F-^ from lattice simulations. To visualize the effect of partially 
quenching we compare the partially quenched two point function to that for Nf = 2 in figure [2j 
When the two masses are identical m„ = rud = m the expression for the partially quenched two point 
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Figure 2: The partially quenched two point function (|3.49p compared to the unquenched two point 
function (j3.36p . Here plotted as a function of the Di eigenvalue keeping the D2 eigenvalue fixed at 
y = 4 and setting m„ = 3 and rhd = 5. The left hand figure is for 5 = 0.4 and the right hand figure is 
for 5 = 1.0. Notice the change of scale on the vertical axis. 



function becomes 



(2) conn r ^ ^\ » ^ 



l+{y,x) Jy{x) xJy+i{x) 



-I^{x,irh) Jy{irh) irhJu-\-i{irh) 
-'jy{x,irh)' Jy{im)' (imJu+iiirh))' 
i~{x,y) e'^'/^My) e-''/^Gu{y,6) 



Jy{irh) imJu+i{im) 
Ju-i{im) imJu{irh) 



(3.50) 

where the prime indicates the derivative with respect to im. 

The two remaining building blocks are given as follows. From the definition (I2.2ip we have from 
eq. SJMh 



h'^Ii{xi,X2) 



det 

hN+iixj + ml){xl + ml) 


HN+i{xi,imu) PN+ii-ml) PN+2{-ml) 
HN+i{xi,imd) PN+ii-ml) PN+2{-ml) 

HN+lixi,X2) PN+lixl) PN+2ixl) 


h^^f\xl+ml){xl+ml) 


PNi-ml) PN+i{-ml) 
PNi-ml) PN+i{-ml) 





(3.51) 



The microscopic limit is readily obtained: 



lim h''^\^{xi,X2) 



det 

2Nxl^^{xl + ml){xl + ml) 


T°(xi,imu) Juiirhu) irhuJu+iiirhu) 
I°{xi,irhd) Ju{irhd) irhdJu+i{irhd) 

I^{xi,X2) Jv{x2) X2Jv+l{x2) 


x^{xl + rhl){xl + ml) 


det 


Juiirhu) irhuJu+iiirhu) 
Ju{irhd) irhdJu+iiirhd) 





(3.52) 
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We note that it is /i/-independent. From eq. (I3.43P and the definition (I2.22P we have 



det 



-ffj^ii (2/1,^2) 



hN+i 



KN+iiyi,imu) PN+i{-ml) PN+2{-ml) 
KN+i{yi,imd) PN+ii-ml) PAr+2(-m^) 

HN+i{yi,y2) XN+i{y2) XN+2{y2) 



leading to 



hm HP,,{yi,y2) = 2N^ 



N^oo 





det 


PN{-ml) PN+i{-ml) 
_ PNi-mj) Pjv+i(-m^) _ 




det 


X+(yi,im„) Ju{irhu) irhuJu+iiirhu) 
I^{yi,imd) Juiiind) irhdJu+i{irhd) 


¥1 


det 


Jy{irhu) irhuJu+iiirhu) 
Ju{irhd) irhdJu+iiirhd) 





(3.53) 



(3.54) 



The two densities are given by 



^(o!i)(y) 



(3.55) 
(3.56) 



They differ because of the asymmetric flavor content, the x eigenvalues not coupling to fl2- Eq- (|3.55p 
agrees with the known result [9], as well as eq. (j3.56p when setting /i2 = 0. 

We are now ready to state the most general correlation function for Ni = 2 and = flavors. 
Upon using the above results, taking out common factors of the determinant we obtain in terms of 
the abbreviations eq. p.l2p 



T°(xii,imrf) 

I (Xj-j^ , Xi2 ) 

I^{yj^,irhu) 
I+{yj^,imd) 

{yji 1 ) 



= n ^» n yi 

i=i j=i 

Ju{irhu) irhuJu+iiirhu 
Ju{irhd) irhdJv+i{irhd 



Juiirrid] 



Juiirhd) 



irhuJu+iiirhu] 
irhdJu+iiirhd) 

{,Xi2 ) 



irhuJv+iiirhu 
irhdJv+iiirhd 

iy{xi^,imd) 
—I [xi-^,yj2) 

I+{yj^,imu) 
I+{yj^,irhd) 

^iyji,yj2) ' 



-| ~n—k 



Juiirhu) 
Juiirhd) 

6 ^ Ju {yj2 ) 

Juiirhu) 
Juiirhd) 

•^-'^I'Uyn) 



irhuJu+iiirhu) 
iTJidJu+iiirhd) 
e-^y'Guiy,2,S) 

irhuJu+iiirhu) 
irhdJu+iiirhd) 

e-^'I'^GuiynM 



(3.57) 



As expected we verify that all correlation functions of eigenvalues Xj are /t-independent and agree 
with the results of [9]. 



4 The non-chiral Two-Matrix Problem 

Having discussed the chiral Unitary Two-Matrix Theory in great detail it is a small step to go through 
the same analysis in the ordinary Unitary Two-Matrix Theory. As mentioned briefly in the introduc- 
tion, this ensemble is relevant for QCD in three space-time dimensions, and may have applications in 
condensed matter physics as well ^16j . 
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It is straightforward to consider two arbitrary chemical potentials m and ^2, as in the chiral 
case treated above. However, for simplicity we restrict ourselves to just imaginary isospin chemical 
potential, i.e. the case where /xi = — /i2 ^ We present the most general case with an even number 
of flavors, 2Nf = 2A^-|_ + 2N^. The techniques used here can also be applied to the case of an odd 
number of flavors. The masses are pairwise equal in magnitude but of opposite sign [14:\ IT5] 



m+i, . . . ,?Ti+Ar^, -m+i, . . . , -m+AT^, and . . . , m_Ar_ , -m_i, , 



-m-N- 



(4.1) 



Both positive and negative mass fermions are coupled to imaginary isospin potential. In intermediate 
steps also an odd number of flavors will appear. The Dirac operator then remains anti-hermitian as 
in four dimensions, and the corresponding Random Two-Matrix Theory reads 



^(2A^/)({m+},{m_}) 



2N+ 2Af_ 

dM^exp [-NTr (^^ + ^2^] det[P+(/i) + m+/] [| det[P_(/i) + m_g](4.2) 

/=i 9=1 



where both <I> and ^ now are N x N hermitian matrices, and 



(4.3) 



A related non-hermitian matrix model aimed at real baryon chemical potential was considered in ref. 
|32j . Here we again restrict ourselves to the Gaussian case for convenience. The same remarks we made 
about non-Gaussian potentials from the previous section also apply here. For the Riemann-Hilbert 
problem to be solved in order to prove universality in the quenched case we refer to )33) . 
We next perform the same change of variables as in eq. ()2.3p to obtain 



27V4 



2Ar_ 



Z^'^^f\{m+}, {m_}) ~ jd^+d^- \[ det[i$+ + m+/] W det[i$_ + rri-g] exp 
with a potential 



9=1 



-^Try(cl>+,<D_) 



(4.4) 



+ 1 /^2 , ^2 



( 



$1 $1 + 2 



1 



-<i) , <i)_ 



(4.5) 



Going to an eigenvalue representation can be done exactly as in refs. \2b\ I26j. Because of the pairing 
in positive and negative masses the integration measure takes the following form 



.00 ^ / ^+ ^- 

z(2^/)({m+},{m_}) = / X{\dxdy^Xl{xj + mlf)X{{y} + 

•'-°°i=i\ f=i 5=1 

,2 



m_ 



xA^({x})A7v({y})det 



(4.6) 



Due to the presence of determinant factors in the measure this is a generalization of the ordinary 
(quenched) Gaussian Unitary Random Two-Matrix Model considered in the literature [25 1 [26 | [T3 l [29 ] . 
We seek biorthogonal polynomials satisfying 



dxdy w('^f\x,y) pi^^^V) Qf'^'hy) = hS^'^hnk , 



with respect to the weight function 



N- 

11(2/' + 

3=1 



mig) e e ' 



(4.7) 



(4.8) 
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The only difference to the previous section is that the /-Bessel function is replaced by an exponential 
and that we integrate over the full real line. Our polynomials are now polynomials in x and y, and 
not in their squares. Furthermore all factors (x^ + are counted as two flavors with masses ^imj, 
instead of one in the previous section. 

Taking into account these differences all definitions of the transforms eq. (j2.19p and the four 
kernels eqs. (j2.20p - (|2.23p carry trough and we shall not repeat them here. The eigenvalue correlation 
functions defined as 



iV!2 



N N 



{N - ny.{N - k)\ Z(2^/)({m+}, {m_}) 



°° j=n+l l=k+l 



i=l \f=l 3=1 / 

(4.9) 

are given by the same formula eq. (|2.25p in terms of the corresponding non-chiral kernels. 

4.1 The quenched case 

The quenched case corresponds to omitting the explicit factors of determinants from the measure, 
while still performing the change of variables (j4.3p . We then get the standard Random Two-Matrix 
Model, with a coupling between the two matrices ^± as given by eq. ()4.5p . We will need these results 
later when expressing the unquenched results in terms of the quenched ones. 
The biorthogonality condition (14. 7[) becomes 

/ dxdy e ^ ^P„(a;)P^(y) = . (4.10) 

Because of the weight being symmetric in x and y we have that Qn{y) = Pn{y)- It was noted indirectly 
in refs. [34j and [28j that this condition is related to Hermite polynomials, although precise details 
were not given there. This can also be seen along the same lines as in appendix |Al In our conventions 
we find for the monic polynomials 



p„w = 2-"(i5i±i)'/f„y^.). (4.11) 

For the normalization constants we obtain 

K = f^rrrr n\ {I - ^l^T . (4.12) 

The factorization of the weight observed in Appendix[B]also applies here. The weight w{x, y) is given 
by the sum over the biorthogonal polynomials 



w{x, y) = exp 



+ 1/2 2x 7*^1 - /^^ 

-iV^(x^ + y^)+iV^x, 



(4.13) 



using the Mehler formula for Hermite polynomials. The individual factors are obviously wi{x) 



W2{x) = e '^^+1 
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The functions Xn{x) can be found explicitly by integrating the polynomials. Alternatively, we 
could also use the results from the previous section by setting u = ±^ and using the relation (|4.18p 
below. The /-Bessel function then becomes a hyperbolic sine or cosine and the two terms add up to 
an integral over the full real line. The result is 



= e-^-H., .l^ . . (4.i4) 

again being proportional to the polynomials. The three different quenched kernels thus read as follows: 

= -fff/fo,!) (4.16) 
^ g^il^..(/=^.)..(/=^.)e-*--,4.., 

With this result we can now obtain all spectral correlation functions via eq. (j2.25p . 

To take the scaling limit, we make use of the relations between Hermite and Laguerre polynomials, 

H2n{x) = (-l)"22"n! 
H2n+i{x) = (-l)"22"+in!x4i/2)(x2) , (4.18) 



and the asymptotic behavior of Laguerre polynomials (j3.6p . This gives 

H2n{X) (-l)"22"n! n-^^cos(V4^X) 

H2n+i{X) ^ (-ir22"+in! ^sin(V4^X) . (4.19) 



If we rescale ^/2Nx = x, V2Ny = y leading to the customary level spacing ^ we arrive at the following 
asymptotic for the kernels: 

lim Hn (y ^ 



lim Mat ( X = = ^H— ) = 2/1^/^ I' {x,y) , (4.22) 



where we have introcued the following notation 



1 ri 2 2 

I^{x,y) = - ds e^'^^''\os{s{x - y)) 
71" Jo 



l\x,y) ^ - I ds cos{six-y)) = Ifl^^i^^l . (4.23) 

TT Jo TT [X -y) 
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The arguments can be real or imaginary. In the asymptotic we have replaced the sums by integrals 
and used again eq. (|3.8p with t = k/N. For convenience we have kept the scaling 2Nfi'^ = p? as in 
the previous section. The weight function has the following limit: 



lim w I X 



X 



V2N y/2N 



exp 



{x - yf 



We will also use the following abbreviation in this section: 



i {x,y) 



1 



exp 



1 

vr Jo 



ds e ^ cos(s(x — y)) 



(4.24) 



(4.25) 



Putting all together we obtain from the definition (j3.15p with (2A^) " ^ replaced by (V^N) " 

^^{yji^Xi^) I°{yj^,yj2) 



P{n,k){{x},{y}) = det 



(4.26) 



after taking out common factors from the determinant. As a check we obtain that the following 
densities are constant and //-independent: 



P{i,o)ix) = P{o,i)ix) 



1 

vr 



(4.27) 



and we also recover the known correlations of the one-matrix model in terms of the sine-kernel for 
P{nfi){x) and P(o,k){y) respectively. 

4.2 Two AND MORE FLAVORS 

In this subsection we outline how the insertion of flavors can be reexpressed in terms of quenched 
quantities as in the last section, and we will also give some examples. 

The biorthogonal polynomials enjoy again a determinental form as in eq. (|3.21|) 

P^'^'^hx) = (-i)'=(det[P+(M)-ix])^(,^,) 



and similarly for 



xAfc({x})Afc({y})det 



N+ N- 
Yl ( dxjdyj{xj - x) Ylix'j + m\f) J|(y| + 

/ 9 
I-/ 

2\i 



exp(jV ^^ .^ Xiy^) 



Qk "(y) = (-i)"(detP_(^)-iy)^(.^,) 
The kernel of both polynomials is given by 

■(det[P+(/^) - ix] det[P_(/i) - iy])^(2Nf) 



(4.28) 
(4.29) 



K^N^(\y,^) 



h 



{2Nf) 
N 



1 



"T / n dxidyi{xi - x){yi - y) W{xj + m\f) W{yf + mi 

f) J-co ■ \ f n 



xA^({x})A7v(M)det 



exp(iV 



2/x2 



'Xiyj J 



f 9 



(4.30) 



22 



with the norms given in eq. (j3.25p . We recah that when applying appendix [O ah different masses in 
eq. (j4.ip are counted as individual ffavors, positive and negative, leading to one column or row in eq. 
(1C.5|) each. The remaining kernels can be obtained from the definitions eqs. ()2.2ip - (j2.23p . The most 
general correlation function then follows from eq. (j2.25p . 

All these quantities, the transformed polynomials and the remaining kernels can be in principle 
be expressed in terms of partition functions. They are given by eq. (|C.5p : eqs. (j3.23p . (j3.26p . and by 
inserting the kernel, (j3.26p . in the definitions (j2.2ip - (j2.23p . with the obvious modification of dropping 
the V terms and writing polynomials with arguments of power one. However, there is a subtlety here 
as the partition function in eq. ()4.28p contains an odd number of flavors. There are two such partition 
functions available I14[ [TSl 



{2Nf+l) 



U(2Nf+l) 



z. 



odd 



U{2Nf+l) 



) 

df/sinh[ySTr(7W[/r5C/t)] , 



(4.31) 
(4.32) 



with Fs =diag(lAr^,lAr^+i), being even and odd in the masses, respectively. Consequently the even 
polynomials P2k are given by eq. (j4.3ip whereas the odd ones -P2fc+i are given by eq. (j4.32p . Similarly 
special care has to be taken when adding single unpaired flavors to each of the sectors A^-|_ and as 
in the kernel eq. (j4.30p . and we refer to p!5] for the relation between partition functions of chiral and 
non-chiral theories at /U = 0. Although similar massive partition functions with chemical potential 
have been computed from complex matrix models [32j and effective field theory [21] we do not further 
pursue this route. Alternatively we can explicitly compute the large-A^ limit from our biorthogonal 
polynomials and kernels, as we will show in the following example. 

Let us consider the simplest case with 2N^ = 2 flavors of masses itm, and no other flavors: 
2A^_ = 0. The corresponding weight function reads 



w 



(2), 



x,y) 



2 2n Ni^xy -Arif£+i(x-2+j;2) 



(4.33) 



The kernel is almost identical to the partially quenched one eq. ()3.43p : 



det 

hN+i 


KN+i{y,im) PN+iiim) PN+2{im) 
Kiy+i{y,-im) PN+i{-im) PN+2{-im) 
KM+i{y,x) Pn+i{x'^) Pn+2{x'^) 


h^J^^\x — im){x + im) 


det 


PN{im) PN+iiim) 
PNi-im) PN+i{-im) 





(4.34) 



The large- limit is easily obtained using the results of the previous subsection, rescaling the mass as 
the eigenvalues y/2Nm = rh: 



lim K 



(2) 
N+l 



{y,x) 



2A^2 


det 

iV2 


I^{y, im) 
T^{y, —irh) 
J0(y,x) 


cosh(m) — isinh(m) 
cosh(m) isinh(m) 
cos(x) — sin(5;) 


(x^ + m?) 




det 


zsinh(m) cosh(m) 
— isinh(m) cosh(m) 





(4.35) 



Here we have chosen to be even, noticing the difference between the asymptotic of Pn and Pn+2 
by a factors of (—1). Odd A^ leads to the same result, as can be seen from interchanging the last two 
columns in numerator and denominator. Similarly the kernel can be seen to be real as it should. The 
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remaining kernels follow with the same modifications from eqs. (I3.45p . (I3.5ip and (j3.53p . Therefore 
we just give the results here: 



det 


I^Oj, irh) 
I^iy, —irh) 
T~{x,y) e 


cosh(?7i) 
cosh(m) 
~^ cos{y) 


—i sinh(m) 
i sinh(m) 
-e~'^^(sin(y) + cos(y)) 


2iV2 


det 


i sinh(m) 
—i sinh(m 


cosh(rn) 
cosh(m) 





(4.36) 



lim h''^^ 



det 

{xl + m^) 


Z^{xi,im) cosh(m) — ?sinh(m) 
Z^{xi,—im) cosh(m) isinh(m) 
I°{xi,X2) cos(x2) -sin(3;2) 


^ det 


i sinh(m) cosh(m) 
— zsinh(?fi) cosh(m) 





(4.37) 



which is /x-independent, and 



V2 N- 



lim H^^\-^{yi,y2) 



We can write the general correlation function as follows: 

.(2) 



det 


x+(yi 

X+(yi, 


im) 
—irh) 
m) e" 


cosh(m) 
cosh(m) 

cos (2/2) - 


—i sinh(m) 
i sinh(m) 
-e-'^''(sin(y2) + cos(y2)) _ 


det 


i sinh(m) 
—i sinh(m) 


cosh(m) 
cosh(m) 





(4.38) 



P\n,k)i{^}n,{y}k) 



det 



I^{xi^,irh) 



-zm 1 



X (xjj , Xi2 ) 

l+{yj^,irh) 
'^^{yji,-^rh) 

iVjl 1 ^12 ) 

From the density 



and from P(o^)(y) in the limit /i ^ we recover the known massive density from the second of ref. 
|14j . Similarly all higher correlators at /i = can be recovered. 

Our second example is for Nj^ = = 1 flavors with masses m+ and m_. The structure is very 
similar to the example in subsection 13.21 but with twice the flavor content: 2Nf = 2 + 2 here. The 
four kernels at finite- A'' are similar to those given by eqs. (j3.3ip . (j3.32p . (j3.37p and (|3.39p . with 2x2 
and 3x3 determinants in denominator and numerator here. For example the first kernel reads 



cosh(m) 
cosh(m) 
cos(xi2) 

cosh(m) 
cosh(m) 
cos(x,o) 



i sinh(m) 
-i sinh(?7i 

-i sinh(rfi' 
i sinh(m) 
-sin(xi2) 



—i sinh(m) 
i sinh(m) 
-sin(xi2) 



cosh(m) 
cosh(m) 

I^{xi^,irh) 
I°{xi^,-irh) 
—X (x , ) 

I+{yj^,im) 
^^{yji,-i'm) 



(4.39) 



cosh(m) 
cosh(m) 
"'^ cos(y 

cosh(m) 
cosh(m) 
"'^^ cos{yj. 



32/ 



n) 



—i sinh(rfi) 
i sinh(m) 
-'^'(sin(yj2) + 2/i2 cos(yj2)) 



—i sinh(m) 
i sinh(m) 
(sin(yj2) + 2/i2 cosii/j^)) 



lim 

Af^oo 



1 o-(2) 

V2N ""^^ 



(x, x) 



(4.40) 
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K 



(2+2) 
N+l 



1 



det 



X det 



with an asymptotic limit 



Kiy^i{—im^,im^) Kiy^i{—im-,—im-^-) 

KN^i{—im^,im+) Kj\f^i{—im^,—im^) Kj\f^i{—im^, x) , (4.41) 
KN+i{y,im+) KN+i{y,-im+) KN+i{y,x) 



lim 

N^oo 



2N' 



1 



X det 



det 



Z^{—irh-,im^) I~^{—im-,—irh^) 

[irh-, —irh^) I^{im-,x) 
I'^(—im-,im+) {—irh- , —irh+) Z~^{—irh-,x) 
Z+{y,im+) Z+{y,-im+) Z+{y,x) 



(4.42) 



Tlie remaining kernels can be obtained similary, and we only give their asymptotic expressions: 



hm M 



(2+2) 
Af+1 



{x,y) 



{x^ + rh^){y^ + mi) 
2iV2 



2A/7r/i- 



det 



X det 



hm hs:?{.ux,) 



(xf + m ^ 
(f 2 + 1TI- 



Z^{im-,iih-^-) Z^ (irh-, —1771-^) 
X+(— im_,im+) Z^{—im-,—irh^) 

Z'^{irh-,irh^) Z'^{im-,—irh^) Z^{im-,y) 

Z^ {—im-,irhj^) Z^{—im-,—imj^) Z^{—im^,y) , (4.43) 

X''(x,im_|_) Z^{x,—im^) Z^{x,y) 

2 



-V2 N- 



1 



det 



X det 



Z~^{im-,im^) X^(im_, — im_|_) 
Z^^—irh-jirh^) Z^{—irh-,—im^) 

Z^{irh-,irhj^) X+(im_,— im+) X^(im_,X2) 
Z^ {—irh^,irhj^) X+(— zm_, — im_|-) X^(— irfi_, X2) 
Z^{xi,irhj^) Z^{xi,—im^) Z^{xi,X2) 



(4.44) 



lim H 

N^oo 



(2+2) 
N+l 



(yi,y2) 



(y^ + ml) 



{yf + in 



V2 N- 



1 



det 



X det 



Z^(im-,im-^-) Z^{im-, —irh^) 
Z~^(—irh^,iih-^-) X+(— im_, — frTi+) 

Z^{im^,imj^) Z^{irh-, —irh-^-) Z^{irh^,y2) 
Z^ {—irh-,i'm+) Z^{—im-,—im+) Z^{—irh-,y2) 
Z+{yi,im+) Z+{yi,-im+) X°(yi,y2) 



(4.45) 



(2-|-2) 

The general correlation function ({x}„, {yj/t) can be easily obtained from the above. As a check 

we have again compared the density 



lim 



(4.46) 



in the limit /i — > to the known result [14j, finding agreement for nondegenerate m+ and m_ 
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5 Summary 



We have introduced a chiral Two-Matrix Gaussian Unitary Ensemble in Random Matrix Theory, and 
solved it explicitly both at finite matrix size and in the — > oo scaling limit. The scaling theory 
describes the Dirac eigenvalue correlations of QCD in the e-regime when subjected to imaginary 
isospin chemical potential. In QCD language the first non-trivial mixed correlation function has a 
strong dependence on the pion decay constant Ft^. This makes it possible to use such correlation 
functions to extract from lattice gauge theory simulations [BE]- Previously two special cases had 
been considered from the viewpoint of effective field theory. Here using the Random Matrix Theory 
formulation we have computed the general correlation functions. We have confirmed that in the two 
special cases {Nf = 0, 2) our general expressions from the present Random Matrix Theory calculation 
agree with the computations from effective field theory [H [2] . 

A particularly useful by-product of the present work is the "partially quenched" case, where the 
Dirac operators that define the theory do no couple to chemical potentials. One can still measure 
the eigenvalue correlations of the Dirac operators Di and D2 on such ensembles, and they will again 
show a very non-trivial dependence on the pion decay constant F-,^. From the point of view of a 
practical implementation in lattice gauge theory simulations this is the ideal set-up. The corresponding 
calculation in the effective field theory seems technically very difficult without some new short-cut or 
trick. Here the Random Matrix Theory provides a novel result that is ready to be used to determine 
the pion decay constant i^V from lattice simulations. 

We have also explicitly found the finite- A'^ solution for the ordinary Two-Matrix Unitary Ensemble 
including mass terms. This extends the existing literature. Our results for the corresponding scaling 
limit are also new. These results are of relevance for QCD in 3 space-time dimensions. We have also 
shown directly how to recover the ordinary One-Matrix results from the Two-Matrix theory. 

Acknowledgement: We thank U. Heller, E. Kanzieper, B. Svetitsky, D. Toublan and J. Ver- 
baarschot for discussions, and we in particular thank B. Svetitsky for suggesting a computation of the 
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DE-FC02-01ER41180. 

A Biorthogonal Laguerre polynomials 

We wish to show here that for Nj = the integral representation eq. (j3.2ip immediately leads to the 
fact that the quenched biorthogonal polynomials are given by Laguerre polynomials as stated in eq. 
(j3.ip . To that aim we will integrate out one of the matrices being Gaussian to map the polynomial 
to the standard matrix representation of the orthogonal Laguerre polynomials in the Gaussian chiral 
one-matrix theory. There the following fact is known: 




det 



ix i<I> 
z<I>^ ix 



exp [-AfTr^>t$l 



(A.l) 
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leading to 

Pk{x^) ~ Li'\Nx') (A.2) 

as orthogonal polynomials with respect to the weight w{x) = x^'^'^^ exp[— A^x^]. Here we have ommited 
all normalisation factors. The matrix $ is complex and of size k x {k + v). 

This can be compared to the representation of our biorthogonal polynomials eq. (|3.2ip : 

Pk{x^) = (-)^(ix)-^^ /"d$id$2det(Pi +ia;)exp[-W($i,«>2)] (A.3) 

with the potential given in eq. ()2.8p . Since Pfc(x^) is given by the expectation value of a determinant 
depending on only one of the matrices, <&i, the matrix integration over $2 decouples. Completing the 
square and integrating out ^2 we obtain an expectation value only with respect to ^>i, with potential 

V{^i) = (ci - d^/c2) Tr$l«>i = Tr$l$i . (A.4) 

Thus we have arrived exactly at eq. (lA.ip with N + /if). Consequently our biorthogonal 

polynomials are Laguerre polynomials as was claimed in eq. (|3.ip . The same argument immediately 
leads to Qk{v^) being Laguerre polynomials as well, integraing out 

It is interesting to compare this to the orthogonality of the Laguerre polynomials in the complex 
plane which arises in the solution of the chiral two- matrix theory for a real chemical potential [lOj . 
A proof of the orthogonality was given in [23] by direct computation. An alternative proof can 
be formulated since the orthogonal polynomials in the complex plane also enjoy a determinental 
representation. Even though the determinant contains both matrices of the two-matrix theory one 
can still find a way to shift the variables of integration such that one of the matrices can be easily 
integrated out leaving a form that again leads to Laguerre polynomials [35]. However the expansion 
method (see Appendix [B]) can not be directly applied since the biorthogonal polynomials Pk{z) and 
p*k{z) are not functions of independent variables. 

B Biorthogonal polynomials from the expansion of the weight 

Consider two sets of polynomials, Pk{x) and qk{y), orthogonal with respect to the two weights wi(x) 
and W2{y) respectively 

J dx wi{x) pk{x) pe{x) = fk Sk/ 

j dy W2{y) Qkiy) qeiy) = Qk h,t ■ (B.l) 

Then it follows automatically that Pk{x) and qk{y) are 6i-orthogonal with respect to the weight w{x, y) 
defined as follows 

00 ^ 

w{x,y) = wi{x) W2{y) 7 Pk{x) qk{y) ■ (B.2) 

The corresponding biorthogonality relation is 

J dxdy w{x,y) pkix) qe{y) = hk 6k/ ■ (B.3) 

While not all weights can be written in the form ()B.2p we will show below that the quenched weight 
(given in (j2.16p with Nj = 0) in fact can. 
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Using these conventions one can easily obtain the transformed functions 

Xkiy) = J dx w{x,y) pk{x) = (hk/gk) W2{y) qk{y) 

Xk{x) = J dy w{x,y) qk{y) = {hk/ fk) wi{x) pk{x) (B.4) 



and ah the kernels 



N-l ^ 

KN{y,x) = IT 1k{y) Pk{x) 

N-l ^ 

Hn{xi,X2) = Wi{xi) ^ — Pk{xi) Pk{x2) 

k=0 

N-l ^ 

HN{yi,y2) = W2{y2) ~ ik{yi) qk{y2) 



k=o 9k 



N-l 



MN{x,y) = wi{x) W2{y) ^ Pk{x) qk{y) ■ (B.5) 

fc=o J 

Interestingly we see that the weight function has the same form as the M kernel just with different 

limits in the sum. We can then write one of the terms appearing in generalized "Dyson Theorem" 
([2:25]) as (for Nf = 0) 

00 ^ 

w{x,y) - MM{x,y) = wi{x) W2{y) ^ ^ Pk{x) qk{y) ■ (B.6) 

k=]\l Jk9k 

In order to apply this to the current problem we make use of the known expansion of in terms 
of Laguerre polynomials giving (for Nf = 0) 

u;(o)(x,y) = {xyy+^ e-^^'''''+'^y'U^{2dNxy) (B.7) 

= {NdYr'^+^xyf^+^e-^-^^^^'+^^y''^ £ ""-^ ~ Y, L'^^KNtcix^) U:\NTC2y'') 

„=o \^ + ^)- 

with r = 1 — dp /{ciC2)- Defining the monic polynomials 

Pk{x) = k\{-NTCi)-^L^^\NTCix^) 

qk{y) = k\{-NTC2)-''Li'\NTC2y^) (B.8) 

and the weights 

wi{x) = 3;2-+ie-^-i-' 

W2{y) = y2i.+lg-JVrC2S;2 ^^g^ 

with support on the positive real line we get the normalisation of the polynomials 

k\{k + iyy. 

~ 2(iVrci)2^'+'^+i 
k\{k + vy 
2{Ntc2) 

From eq. ()B.7p we can then read off the norm of the biorthogonal polynomials 

_ k\{k + v)\{Ndf^+'^ 

The same considerations also apply to the non chiral case in terms of hermite polynomials, cf. eq. 
(I4l3l) . 



9k = ^u - (B-io) 



28 



C All partition functions 

In this appendix we derive the most general partition function, containing a different number of A^i 
flavors with masses rufi and flavors with masses mf2, where A'^^ = iVi + This result is 

particularly useful as both the unquenched polynomials Pjf^^^ and Q^f^^^ as well as their kernel can 
be expressed in terms of these partition functions, see eqs. (|3.23p and (|3.26p . 

Following the same steps as in section [2] the partition function has the following eigenvalue repre- 
sentation: 

■^0 i=\ \ f=i f=i 

X A^({x2})A^({2/2})det[/,(2d7Vx,y,)]e~'^^»^'^"^+'^'^'^ ■ (CI) 

Because all eigenvalues are integrated out the symmetry of the Vandermonde determinants can be 
used to bring all A^! terms in the sum over permutations from the determinant of the /-Bessel function 
into diagonal forrr0 



^r^\Wi},{m2}) = M{{m,,2})N\ ri[[dx,dy,ixmr'''t[{xj + m%)l[{yf + mj. 

° i=i V f=i f=i 

N 

X A^({x2})A^({y2})]^4(2diVx,y,)e-^S.(^^^?+^^^?) . (C.2) 

i=l 

Thus we can now interpret the /-Bessel function as being part of the weight. In this form we can 
apply a theorem that was proven in [20] for matrix theories with complex eigenvalues, expressing such 
partition functions (also called characteristic polynomials) in terms of the quenched polynomials. We 
can apply this result by formally identifying z ^ x and z* ^ y, thus relating our polynomials to the 
polynomials V there as follows: 'Pfc(-z^) Pjfc(x^) and Vk{z'^)* Qkiv^)- Note that unlike in [20| . 
where the two sets of polynomials are related by complex conjugation, our polynomials P and Q are 
no longer related by a symmetry in general. However, this property was not used in the proof of the 
theorem we wish to apply. 

The quenched kernel made of the monic polynomials without the weight, 

(C.3) 

is obviously self-contractive, due to the biorthogonality: 

dxdy w'^°\x,y) KN{y,u)KNiv,x) = Kn{v,u) . (C.4) 



This makes Dyson's Theorem for integrating over determinants of such kernels applicable, which was 
used in the derivation in [20]. We therefore arrive at 



X det |^KAr+Ar2(im/2,imfci) PN+N2{{'i'mkif) . . . PN+Ni-iHimki f) (C.5) 



^This is no longer true in correlation functions iin.fc when less than A*' of each eigenvalues are integrated out. 
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where we have given the result for A'^i > N2. The first N2 cohimns are given by kernels and the last 
Ni — N2 columns by monic polynomials of increasing degree. The corresponding result for N2 > Ni 
then obviously follows by exchanging A^i <-> N2 and irifi -(r^ mf2 ib. the above formula. 

The fact that imaginary arguments appear is related to the fact that mass terms are positive, 
ni(a:;i + ml). Polynomials and kernels of real arguments would appear, if we had considered charac- 
teristic polynomials, such as Jlil^? ~ ^^)- 
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